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ABSTRACT 

Aims. The aim of this work is to understand to what extend gravitational interactions between the stars in high-density young 
stellar clusters, like the Orion Nebula Cluster (ONC), change the angular momentum in their protoplanetary discs. 
Methods. Two types of simulations were combined — N-body simulations of the dynamics of the stars in the ONC, and angular 
momentum loss results from simulations of star-disc encounters. 

Results. It is shown that in a star-disc encounter the angular momentum loss is usually larger than the mass loss, so that the 
disc remnant has a lower specific angular momentum. Assuming an age of 1-2 Myr for the ONC, the disc angular momentum in 
the higher density region of the Trapezium is reduced by 15-20% on average. Encounters therefore play an important part in the 
angular momentum transport in these central regions but are not the dominant process. More importantly, even in the outer 
cluster regions the angular momentum loss is on average 3-5%. Here it is shown that a 3-5% loss in angular momentum might be 
enough to trigger gravitational instabilities even in low-mass discs - a possible prerequisite for the formation of planetary systems. 
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1. Introduction 

Discs are seen as a prerequisite for the formation of plan- 
etary systems. There is increasing observational evidence 
that most, if not all, stars are initially surrounded by discs. 
For example, Lada et al. (2000) found that 80-85% of all 
stars in the Trapezium are surrounded by discs. 

One visually compelling feature in the numerical study 
of these discs is the development of spiral arms either in 
self- gravitating high -mass discs or as an effect of encoun- 
ters lPMzHer20o3) . In a few cases where discs have been 
observed with a sufficiently high spatial resolution, the 
existence of suc h spiral structures ha. v e been observation- 
ally con firmed IClampin et al. 2003t ICorder et al. 20051 
iLin et al . 2006, for example). Spiral arms are always con- 
nected to mass and angular momentum transport. Our 
study of the importance of the encounter-induced angular 
momentum transport in an cluster environment has two 
motivations: angular momentum transport (a) in the late 
stages of accretion in star formation and (b)in spiral arm 
formation as prerequisite for gravitational instabilities and 
fragmentation in planet formation scenarios. 

The former is connected to the problem that the typ- 
ical o bserved angular mo mentum of the cloud cores is so 
large ijBodenheimer IQQSiD that the star would be spin- 
ning at velocities higher than its break-up velocity if no 
means for angular momentum loss would be present. One 



way out is angular momentum transfer to the disc and 
then from the inner disc regions outwards. Many differ- 
ent mechani sms have been proposed for the latter (see 
iLarson 20021 and references therein). Although often dis- 
carded as probably not important, encounter- induced an- 
gular momentum transport is still often cited as possible 
mechanism. To clarify this point we study disc angular 
momentum loss in a cluster environment for the first time 
in a quantitative, systematic manner. 

The second aim of this investigation is to achieve a 
better understanding of how the cluster environment in- 
fluences the development of spiral structures in the discs 
through encounter-triggered angular momentum trans- 
port. This might be of special importance as thus such 
spiral structures are a prerequisite in the theory of planet 
formation through gravitational instabilities. 

The Orion Nebula Cluster (ONC) was chosen as the 
model cluster for these investigations because it is one of 
the best-studied regions in our galaxy, and observational 
constraints significantly reduce the modelling parameters. 
In addition, its high density suggests that stellar encoun- 
ters might be relevant for the evolution of circumstellar 
discs. This is supported by recent investigations by Olczak 
et al. (2006) indicating that the effect of encounters has - 
at least for the disc mass loss - previously been underes- 
timated. 
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The investigation involves three steps — i) the deter- 
mination of the encounter parameters in the ONC dynam- 
ics, ii) the angular momentum loss in isolated encounters, 
and iii) the combination of both results to investigate the 
angular momentum loss of the discs in the cluster. 

For the first step we use the cluster simulations de- 
scribed in Olczak et al. (2006) and star-disc simulations 
in Pfalzner et al. (2006) to determine a parameter list for 
the encounters for a model ONC in equilibrium and ex- 
pansion, i.e. virial ratios of Qyir=0.5 and Qyir—1. The 
results are summarized in Section 2. 

This is followed by an investigation of the angu- 
lar momentum loss in star-disc encounters in Section 3. 
There have been earlier investigations which looked at 
angular momentum transport in isola t ed star-disc en- 
counters ("Ostri ker IQol iLarwood 1997t iHall et al. 19961: 
rPfalzner 2004} predominantly in prograde, coplanar en- 
counters for low-mass discs covering only a small param- 
eter range. Under these limitations simple N-body simu- 
lations suffice, and hydrodynamic effects and self-gravity 
within the disc can be largely neglected. Following this ap- 
proach one has to keep in mind that the results have to be 
considered as upper limits (see discussion in Olczak et al. 
(2006) for the disc mass loss). We extended the parame- 
ter range covered in Pfalzner et al. (2005a) to larger mass 
ratios between the interacting stars and penetrating en- 
counters to cover all situations necessary when modelling 
the ONC. In this work we neglect photo-evaporation as 
disc dispersal mechanism (Scallv & Clarke 2001), but this 
should be included in a future investigation. 

In Section 4 the results of Section 2 and 3 are combined 
to determine the disc angular momentum as a function of 
time. In Section 5 it is shown how the angular momentum 
mass loss is influenced by the assumptions made and a 
discussion is given of how this angular momentum loss 
can trigger gravitational instabilities. 

2. Model and Cluster simulations 

We combine simulations of the dynamics of the ONC 
determining the interaction paramete rs of clos e encoun- 
ters between stars in the cluster (Olczak et al. 20 06') and 
with result s of isolated star-disc encounter simulations de- 
scribed in jlPfalz ner et al. 2006|) . In contrast to these pre- 
vious studies, the emphasis in the current investigation is 
not on the mass but the angular momentum loss of the 
discs. 

The dynamical model of the ONC presented here con- 
tains only stellar components neglecting gas and the po- 
tential of the background molecular cloud OMC 1. The 
cluster models were set up with a spherical density dis- 
tribution p(r) oc r^^, by randomly generating the masses 
according to the mass function given by Kroupa et al. 
(1993) in a range SOM© > M* > O.OSMq. 0^C Ori was di- 
rectly assigned a mass of 50 Mq and placed at the cluster 
centre. The velocity distribution of the stars was assumed 
to be Maxwellian. 



The O NC cluster sim ulations were performed with 
NBODY6++ ijSpurzem 199^ and modelled for 13 Myr - the 
assumed lifetime of 9^C Ori. The quality of the dynam- 
ical models was judged by comparing to observational 
data at 1-2 Myr, marking the range of the mean ONC 
age. For more details of the selection process, see Olczak 
et al. (2006). The quantities of interest were: number 
of stars, half-mass radius, number densities, velocity dis- 
persion and projected density profile. Here we describe 
mainly the case of the ONC being in virial equilibrium, 
in Section 01 the changes one can expect for an expanding 
cluster will be briefly considered. 

Assuming that only two-body encounters occur and 
that higher-order encounters are negligible, the effect of 
an encounter can be investigated by considering it to be 
isolated from the rest of the cluster. During the course 
of the simulation the information of all perturbing events 
of each stellar disc, i.e. both masses, the relative velocity 
and the eccentricity are recorded. However, the latter is 
not applied as such later on. To limit the parameter space 
only the angular momentum loss of parabolic encounters 
is used. The eccentricity values are utilized to estimate the 
error introduced by such a simplification. 

3. Angular momentum loss in star-disc encounters 

This part of the investigation is mainly based on the ear- 
lier work by Pfalzner et al. (2005a) where a parameter 
study was performed for a star of mass Mi = 1 Mq sur- 
rounded by a disc perturbed by the fly-by of a second star 
of mass When simulating isolated encounter events, 
the observational evidence that most discs in the ONC are 
of low mass m^, i.e. md/M* <C 0.1, reduces the complex- 
ity of the numerical approach in several ways: low-mass 
discs do not significantly influence the encounter orbit and 
self-gravitation can be neglected and the results are scal- 
able to other star masses. The effect of viscosity is in this 
case negligible too, as the viscous timescale is so much 
longer than the timescale of the encounter. Another sim- 
plification is that only coplanar, prograde encounters were 
investigated. According to st udies on incli r ied and ret- 
rograde star-d isc enco unters llOstriker 199^ iHeller 19951: 
iHall et al. 1996: Pfalzner et al. 2005aj) . due to these ap- 
proximations the results can only be interpreted as an up- 
per limit for the angular momentum loss. We will dicuss 
this in more detail at the end of this section. 

The disc is assumed to extend to rd = 100 AU and the 
surface density initially has a l/r-dependence. The angu- 
lar momentum loss induced by the fly-by of a star on such 
a star-disc system depends on various parameters specific 
to the stars and discs involved an d the encounter orbit . 
The above mentioned earlier study l|Pfalzner et al. 2005af) 
was limited to non-penetrating encounters with star mass 
ratios in the interval 0.1 < M2/M* < 5. Here, the pa- 
rameter study is extended to the higher mass ratios of up 
to M^/M^ = 5O.OM0/O.O8M0 which occur in the ONC 
simulations. The additional fractional angular momentum 
loss (AML) values can be found in Table H available with 
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Fig. 1. Interaction range for mass and AML as a function 
of the mass ratio of the interacting stars in units of the 
disc radius. 

the electronic form of this article. The improved fit func- 
tion for (AJ/ J)j„,^; is 
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where = rperi/rdi 



is the periastron separation rperi in 



units of the disc radius rdisc- 

The most important results for the cluster are that 

1. the AML is always larger than the relative mass loss, 

2. the angular momentum of the disc is affected by much 
more distant encounters than the disc mass. 

The latter is illustrated by Fig. ^ where the relative 
distance Vperi/fdisc at which the angular momentum is 
changed by at least 3% (the average error value) is shown 
as a function of the mass ratio M|/M*. Given a suffi- 
ciently large perturber mass, even encounters as distant 
as 20 times the disc radius can reduce the angular mo- 
mentum in the disc by 10% or more without changing the 
disc mass at all. 

A considerable part of the angular momentum is 
carried away with the matter that becomes unbound. 
However, what one is really interested in is to what degree 
the angular momentum per particle is lower in the remain- 
ing disc — only this can eventually facilitate accretion of 
matter onto the star. A measure of this specific angular 
momentum can be defined by 



{^JlJ)raass ^= ( ^ " ^ 



1 - 



Ar 



(2) 



Fig. [2Ji and b show the specific angular momentum ob- 
tained by using Tablenin this paper and Table 3 in Olczak 
et al. (2006) as a function of Mj/Mj* and r^, respectively. 
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Fig. 2. Specific angular momentum after an coplanar, pro- 
grade parabolic encounter compared to that of the initial 
disc shown as a function of a) the perturber mass M| and 
b) the encounter periastron in units of the disc size. 

^ For large mass ratios there exists an lower limit of 0.4 for 
{Jenc/ J)mass Considered periastra. This means that 

in a single encounter, however strong, the specific angular 
momentum can only be reduced to ~40% of that in the 
initial disc. However, such a change of angular momentum 
is achieved over a wide range of interaction parameters, so 
that in repeated encounters the angular momentum of the 
remaining disc could nevertheless be considerably reduced. 

Since the cluster consists of a wide spectrum of star 
masses, these simulation results obtained for Mj* = 1 Mq 
are generalized by scaling the disc radius according to 



rd^.sc = r<j,,,(lM0)y^M*[M0] 



(3) 



In the cluster simulations r^isc (IMq) was assumed to be 
150AU. Scaling the disc size with the star mass seems in- 
tuitively right, but observations are not so unambiguous: 
Vicente & Alves (2005) found a correlation between disc 
diameters and stellar masses using a sample of proplyds 



^ Note, fit formula [H here and Eq. 4 in Olczak et al. (2006) 
are not of high enough accuracy to determine {Jenc/J)^^^^ in 
Eq.|21in the entire parameter range required for the ONC. 
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Fig. 3. a) Percentage of stars in bins with a 5% (lowest 
considered bin) or 95% (highest considered bin) AML in 
the disc as a function of distance to the cluster center at 2 
Myr averaged over 20 simulations, b) shows the same as a 
function of the mass ratio Af|/Af* of the two interacting 
stars. For comparison in both figures the same is shown 
for the mass loss, too. 



replenished with material of very low angular momentum 
so that the specific angular momentum loss (SAML) would 
be underestimated. 

Most encounters in the cluster simulations are not 
parabolic but hyperbolic. In such hyperbolic encounters 
the angular momentum loss is lower because the disturber 
is not long enough in the vicinity of the star-disc system to 
remove angular momentum in an efficient way. However, 
considering only the stars that lose more than 90% of their 
disc mass, the eccentricity distribution has a maximum at 
e « 3. Pfalzner et al. (2005a) showed that for = 1 Mq 
the AML in an e = 3 encounter is about 55% of that of a 
parabolic encounter. 

In a cluster that is not highly flattened it seems rather 
unlikely that the encounter planes are aligned to a high 
degree. Therefore one would expect most encounters to 
be non-coplanar. Pfalzner et al. (2005a) showed that, as 
long as the inclination is not larger than 45° the AML in 
the encounter is only slightly (< 10%) reduced in compar- 
ison to a coplanar encounter. If however the orientation 
is completely random and a 90°-encounter the most likely 
encounter scenario, the mass loss could be significantly 
reduced, a point which would need further investigation. 

It is often argued that consecutive encounters would 
not lead to any significant additional momentum trans- 
port because the outer edge of the disc is truncated by 
the first encounter - the part that is mainly involved in 
the interaction. However, Pfalzner et al. (2004) demon- 
strated that a second encounter can indeed lead to the 
same relative AML than in the first encounter. This 
implies that a succession of distant encounters might 
well be able to transport a higher amount of angu- 
lar momentum outwa rds than previously thought (see 
iMoeckel fc Ballv 20061 too). In future studies, it should 
be investigated whether this holds for the entire parame- 
ter space covered here. 



from Luhman et al. (2000), but not for the data from 
Hillenbrand (1997). However, the present Trapezium is 
probably not in its primordial state as various disc de- 
struction processes have most likely altered the disc sizes 
considerably. 

As mentioned earlier the encounter results summarized 
by Eq. 13 Fig. [5] and Table ^ represent upper limits to 
the angular momentum loss induced by encounters. The 
reason is that they are obtained for prograde, coplanar 
and parabolic encounters only. In addition, it was assumed 
that only one of the stars was surrounded by a disc. The 
latter is somewhat contradicting - assuming that each star 
is initially surrounded by a disc but at the same time con- 
sidering an encounter with a disc-less star. There are two 
reasons why this was done: first, encounters where both 
stars are surrounded by discs are less well investigated and 
second Pfalzner et al. (2005b) showed that the star-disc 
results can be generalized to disc-disc encounters as long 
as there is no mass exchange between the discs. In the 
case of a mass exchange the discs can be to some extent 



4. Angular Momentum Loss in the ONC 

In the following the cluster results are combined with those 
of Section to investigate the encounter-induced angular 
momentum loss of the discs in the ONC. First we deter- 
mine the AML using Eq. 2] The fact that the angular 
momentum loss is usually larger than the mass loss, is 
reflected by the cluster simulations in several ways: As 
Fig. |2t) shows, the number of stars with, for example, 
more than 5% or 95% relative angular momentum loss in 
the disc is somewhat higher than the number of discs with 
a more than 5% or 95% mass loss. Due to the lower stel- 
lar number density in the outer cluster regions both the 
mass loss and the AML decrease with distance from the 
cluster center. However, the AML extends to more distant 
cluster regions — for the 5% case the mass loss extends 
to ~ 0.75 pc , whereas the AML reaches to ^ 1.25 pc. For 
high losses (95%) the mass loss extends only to 0.3-0.4pc, 
but the angular momentum to ^ l.lpc. Fig. |2l3 demon- 
strates that lower relative perturber masses M^/Mj* are 
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Fig. 4. Percentage of stars with a SAML as defined by 
Eq. 121 as function of the stellar mass averaged over 10 
simulations. The bin for 5% and now 55% (as the highest 
obtainable SAML is ^ 60%) loss are shown. 

required to remove at least 5% of the angular momentum 
than mass from the disc. 

The dependence on the distance from the cluster cen- 
ter of the SAML is very similar to that of the AML, 
i.e. whereas beyond 0.3pc no mass loss occurs, AML and 
SAML happen as far out as ^ 1.25 pc. However, some dif- 
ferences nevertheless occur. They are due to the fact most 
massive stars have basically lost all their mass and are 
not considered anymore for this mass-corrected value. In 
addition, the dependence on AI^/Mi is moved to higher 
Af| /Mj*-values (see Fig.QlD). This might indicate that low- 
mass encounters can remove only mass with angular mo- 
mentum, whereas high-mass encounters can do the oppo- 
site. 

Fig. 13 shows that the SAML is reduced most in the 
central region of the cluster (by ^ 15-20%), but through- 
out the cluster there is SAML even if only on a 3-5% per 
cent level in the outer cluster areas. 

As mentioned in Section 3 it is not clear what the pri- 
mordial disc-size distribution in a cluster is. If we assume 
that the disc size does not scale with the disc mass, but as- 
sume initially equal disc sizes (150AU), the SAML is larger 
at all distances (see Fig. |S| . This is most pronounced in 
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Fig. 5. SAML as a function of the radial distance to the 
cluster center averaged over 20 simulations. 
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Fig. 6. SAML for equal and scaled disc size as a function 
of the radial distance to the cluster center. 



the Trapezium (~ 30-40%), but even further outside the 
SAML is increased by ^ 1%. The reason is that the ma- 
jority of stars have masses <1M0. Since an equal disc size 
for all stars is equivalent to a larger disc size for them, this 
results in a higher SAML. 

As the velocity dispersion of the ONC is not well de- 
termined, the next point to investigate is how the SAML 
changes if we consider the cluster to be expanding. Testing 
the Qvir=i case, the angular momentum reductions turns 
out to be larger in the entire ONC with the largest differ- 
ences in the Trapezium region (see Figl?)). 

Fig.|Slshows that the difference in SAML happens early 
on in the cluster development. Here the density in the 
expanding cluster is initially higher in the center leading 
to more and stronger encounters. 
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Fig. 7. SAML as function of radial distance to cluster cen- 
ter. Comparison between cluster in equilibrium and an 
expanding cluster for scaled disc size. 



Fig. 8. SAML as function of time for an expanding cluster 
at O.lpc and l.Opc. 



5. Discussion and Conclusions 

Combining simulations that follow the path of the stars in 
the cluster with angular momentum loss results obtained 
from simulation of isolated encounters, the encounter- 
induced angular momentum loss has been investigated for 
the ONC. The main outcome of these simulations is that 
for most stars in the ONC, encounters reduce the angu- 
lar momentum on average only by 3-5% in the remaining 
disc compared to the initial disc. However, this kind of 
reduction happens even in the outer parts of the ONC, so 
that the disc angular momentum is affected by encoun- 
ters throughout the entire cluster. If one considers only 
the central dense Trapezium region the SAML increases 
to the 10-30% level. This means that although encounters 
in the cluster environment do lead to angular momentum 
loss in the disc, it is not the predominant process for an- 
gular momentum transfer. 

However, it would be premature to declare encoun- 
ters as not important for the disc development. On the 
contrary, encounter-induced angular momentum transport 
may still play an important part in the disc evolution. As 
one can see in Fig.O an SAML in the 3-5% regime, which 
we showed to be present in the entire ONC, is sufhcient to 
trigger pronounced spiral arm structures in the disc. As 
one can see similar structures appear independent of the 
mass of the perturber given the same angular momentum 
change. 

As mentioned above such spiral structures have been 
detected in the few cases where discs have been observed 
with a sufficiently high spatial resolut i on (for exam ple 
IClampin et al. 2003t ICorder et al. 20051 iLin et al. 20 06). 
A fly-by has often been rejected as possible cause of the 
spiral pattern because no stars have been found in the im- 
mediate vicinity of the star-disc system. However, as the 
results in Section 3 show, angular momentum change can 
be induced by relatively distant encounters given that the 



mass ratio between the star and the perturber is large. 
For example, for M2 /Ml—50 an encounter with a perias- 
tron ~14 times the disc radius could create a spiral arm 
pattern similar to those shown in Fig. |51 Considering the 
far range of interaction in many cases, a much larger area 
would have to be considered in a search for possible en- 
counter partners. 

What consequences can such spiral patterns in the disc 
have? The spiral arm pattern above might be the starting 
point for gravitational instabilities (in above simulations 
self-gravity of the disc is not considered, so clumping can 
not occur). So far gravitational instabilities have predom- 
inantly been in yestigated for isolate d high-mass discs (for 
a summary, see lPtirisen et al. 2006|) . In this context grav- 
itational instabilities can produce global spiral arms and 
disc fragmentation into dense clumps and substructures. 
Possibly these dense clumps may become precursors to gas 
giant planets l|Kuiper 195lt ICameron 19781: IBoss 1998|) . 

Modern numerical sirn ulations, beginning with 
l)PaDaloizou fc Savoniie 199lf) . show that non- 
axisymmetric disturbances become unstable for 
Q — Csk/ttGT, < 1.5, where Cs is the sound speed, 
K the epicyclic frequency, G the gravitational constant 
and E the surface density. When the gravitational 
instabilities emerge from the linear regime, they may 
either saturate at nonlin ear amplitude or fragment the 
disc ijPurisen et al. 200^ . What happens depends on 
whether a balance can be reached between heating and 
the loss of disk thermal energy by radiative or convective 
cooling. Although there is agreement on conditions 
for fragmentation, it is an open question whether real 
discs cool fast enough for fragmentation to occur and 
whether these fragments last long enough to contract into 
permanent protoplanets (jPuriscn ct al. 2 006). 

A prerequisite for the gravitational instabilities in iso- 
lated discs is a high enough disc mass that self-gravity can 
carry small perturbations to fragmentation in the non- 
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Fig. 9. Examples of encounters with a (SAML) of 3-5%. 

a) shows an encounter where M2 = 0.1 Mq and Vp = 2.0, 

b) an encounter where M|= 1 Mq and rp = 3.5 and c) 
an encounter where M2 = 5 Mq and Vpd = 5.5. 




Fig. 10. Examples of encounters with a SAML of ~20%. 

a) shows an encounter where M|= 0.1 Mq and Vp = 1.5, 

b) an encounter where = 1 Mq and rp = 2.8 and c) 
an encounter where M2 = 5 Mq and rp = 4.5. 



linear regime. However, observations show that high-mass 
discs seem to be the exception rather than the rule. 

Various mechanisms have been suggested 
for triggering gravitational instabilities, among 
them formation of massive discs from the col- 
lapse of protostars llLaughlin fc Bqdenheimer 1994 
lYorke fc Bodenheimer 1999t IBoss 1998|) . clumpy ki- 
fall onto the disc IjBoss 19981) . cooling from stable to 
unstable state, p erturbation by binary companions 
I Maver et al. 20 0^ and close encounters with star/disc 
systems (,Bof5n et al. 19981 iLin et al 19981). 

Investigating binary systems l|Maver et al. 20051) 
found that the degree of fragmentation depends on the 
mass of the disc and the cooling rate. The simulations 
by Mayer et al.(2005) and Boss (2006) both indicate that 
apparently, in contrast to isolated discs, in binary systems 
lighter discs are more likely to fragment. This counterin- 
tuitive result reflects the fact that intense heating from 
spiral shocks can suppress the formation of permanent 
clumps. A balance between compression/shock heating 
and cooling is required for fragmentation. 

The situation in an encounter is similar to that in 
a binary system, because the secondary star induces a 
non-axisymmetric perturbation onto the star-disc system, 
but with the difference that the perturbation is not re- 
curring. At the moment it is unclear whether encoun- 
ters enhance di sc fragmentation or not. Investigations by 
I Lodatao et all ) have found that discs with long cooling 
timescales are actually rendered more stable by encoun- 



ters. On the other hand, ijShen fc Wadslevf) do get en- 
counter induced fragmentation when they use an isother- 
mal equation of state. In the future it needs to be deter- 
mined by detailed simulations with the right radiation hy- 
drodynamics which amount of angular momentum change 
is required to obtain stable clumps — is it the stronger 
interactions (10-30% angular momentum change) that we 
find in the cluster center or rather the weaker (3-5% an- 
gular momentum change) ones further outside? 

There is actually a second means by which spiral struc- 
tures induced by encounters could support planet for- 
mation. In the classical picture o f planet (for a sum- 
mary see IWeidenschilling fc Cuzzi|) formation sticking of 
dust grains and further agglomeration leads eventually to 
planet formation. Rice et al. 2006 showed that dust can 
accumulated in spiral shocks induced in self-gravitating 
disc. If this happens in the same fashion in encounter in- 
duced spirals, encounters could trigger planet formation 
via dust agglomeration as well. The higher dust density 
and possibly the higher velocity dispersion in the spiral 
arm could assist dust agglomeration. However, the time 
scale of the lifetime of spiral structures is typically only 
1000 to 2000 years which would be very stringent for this 
planet formation scenario. 

From the simulations in this paper we have shown that 
in systems like the ONC such an encounter situation is a 
common event throughout the cluster. Therefore, if plan- 
ets can be formed through gravitational instabilities at all, 
encounters in a cluster environment are a likely trigger for 
their occurance. 
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Table 1. Table of AML 



(AJd/Jd) bound 

for all simulated configurations of parabolic (e = 1) star- 
disc encounters. The first row contains the relative per- 
turbcr masses ili^l/Mj*, the first column contains the rela- 
tive periastra rp/rj- Results from simulations are denoted 
by four digits, the values "0.0" were edited manually. 



